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ABSTRACT 



We present spectroscopic and photometric observations of the optical counterpart to PSR J1911-5958A, a millisecond pulsar located towards 
the globular cluster NGC 6752. We measure radial velocities from the spectra and determine the systemic radial velocity of the binary and 
the radial-velocity amplitude of the white-dwarf orbit. Combined with the pulsar orbit obtained from radio timing, we infer a mass ratio of 
Mpsr/Mwd = 7.36 ± 0.25. The spectrum of the counterpart is that of a hydrogen atmosphere, showing Balmer absorption lines up to H12, 
and we identify the counterpart as a helium-core white dwarf of spectral type DA5. Comparison of the spectra with hydrogen atmosphere 
models yield a temperature T^fi = 10090 ± 150 K and a surface gravity logg = 6.44 ± 0.20 cgs. Using mass-radius relations appropriate for 
low-mass helium-core white dwarfs, we infer the white-dwarf mass Mwd = 0. 18 ± 0.02 Mq and radius /?wd = 0.043 ± 0.009 i?o- Combined 
with the mass ratio, this constrains the pulsar mass to MpsR = 1.40^q Mq. If we instead use the white-dwarf spectrum and the distance of 
NGC 6752 to determine the white-dwarf radius, we find /?wd = 0.058 ± 0.004 A©- For the observed temperature, the mass-radius relations 
predict a white-dwarf mass of Mwd = 0. 175 ± 0.010 Mq, constraining the pulsar mass to MpsR = 1 .34 ± 0.08 Mq. We find that the white-dwarf 
radius determined from the spectrum and the systemic radial velocity of the binary are only marginally consistent with the values that are 
expected if PSR J191 1-5958A is associated with NGC 6752. We discuss possible causes to explain this inconsistency, but conclude that our 
observations do not conclusively confirm nor disprove the assocation of the pulsar binary with the globular cluster 

Key words. Pulsars: individual (PSR J191 1-5958A) - globular clusters: individual (NGC 6752) - stars: neutron - stars: white dwarfs 



1. Introduction 

The equation-of-state of matter at supra-nuclear densities 
together with general relativity implv a m aximum mass 
for a rotating neutron star (e.g. iLattimer & Prakash .2004 ). 
Conversely, a measurement of a high neutron-star mass con- 
strains the equatio n-of-state of this matter at these densities. 
iThorsett & Chakrabartv. (.1999.) found that neutron stars in ra- 
dio pulsars cover only a rather narrow range in mass; 1.35 ± 
0.04 Mq. However, their sample is statistically dominated by 
mildly recycled pulsars in relativistic double neutron-star bina- 
ries. Considerably higher masses (up to ~ 2Mq) are expected 
for millisecond pulsars with low-mass white-dwarf compan- 
ions, since binary evolution predicts that several tenths of solar 
masses of material have been transferred from the progenitor of 
the white dwarf onto the pulsar, spinning it up to the cuiTently 
observed (millisecond) periods (for rev iews, see . VerbunL 1 993: 
IPhinnev & Kulkamill994tLStairsl2004 . 

About 40 of such s vstems are known (see review by 
Ivan Kerkwiik et allbOOSi) . but neutron-star masses have been 
measured for only six of them. PSR J075 1-1-1807 contains 
the heaviest n eutron star known to date and with a mass of 
2.1 + 0.2 Mq jNice et alj. l2005 b) this is the only system for 
which the mass is not consistent with a value near 1 .4 Mq. 



With the exception of PSR J1012+5307, these neutron-star 
masses are determinated from radio timing of the millisecond 
pulsar; either due to the detection of general-relativistic ef- 
fects or due to the detection of secular and annual variations 
because of the motion of the Earth. For PSR J1012-i-5307, 
the neutron-star mass was determined through optical spec- 
troscopy of the white-dwarf companion to the pulsar These 
measurements yield the radial-velocity amplitude of the white- 
dwarf orbit, which combined with the pulsar orbit, determines 
the mass ratio between the white dwarf and the pulsar A 
model-atmosphere fit to the white-dwarf spectrum provides the 
effective temperature and surface gravity of the white dwarf. 
Combining these values with white-dwarf mass-radius rela- 
tions yield the white-dwarf mass and radius and, through the 
mass ratio, the pulsar mass. 

In this paper, we use this method to determine the mass 
of the binary millisecond pulsar PSR J191 1-5958A. This pul- 
sar is in a 20 h, highly circular (e < 10"^) orbit around a low- 
mass companion and located at a projected of fset of 6'.4 from 
the center o f the globular cluster N GC 6752 jP' Amico et all 
1200 iL l2002i) . ID' Amico et"al] 12002) argued that the pulsar bi- 
nary is associated with the globular cluster NGC 6752 (cf. the 
discussion in Appendix IB. 3> . In order to explain the large dis- 
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tance of the pulsar fro m the cluster center (3.3 half-mass radii) 
and the circular orbit, IColni et alJ (l20Q2h investigated several 
possible scenarios. They argued that if PSR J 19 11 -5958 A was 
ejected out of the core of NGC 6752 this may be the result of 
an encounter with a wide binary consisting of two black holes. 

The optical counterpart to PSR J191 1-5958A was discov- 
eredby Bassa etal. (2003b) and confirmed by Ferraro et al. 
( I2OO3I) . It was found that the colours and magnitudes of 
the counterpart are consistent with those of a 0.1 8-0.20 M© 
helium-core white dwarf at the distance of NGC 6752. The 
relative brightness of the counterpart {V = 22.1) and the 
fact that the field surrounding PSR J1911-5958A is not ex- 
tremely crowded, motivated us to obtain phase-resolved spec- 
troscopic observations of the companion of PSR J191 1-5958A 
and determine the mass of the pulsar. In principle, these 
observations can also be used to verify the membership of 
PSR J1911-5958A with NGC 6752 through the systemic ra- 
dial velocity and the white-dwarf radius, which should be con- 
sistent with values expected for a system associated with the 
globular cluster. If the association is confirmed, the accurate 
distance to the globular cluster provides a separate constraint 
on the radius of the white dwarf and thus its mass. 

This paper is structured as follows; in § |3 we describe 
our spectroscopic observations and their reduction, as well as 
the analysis of archival photometric observations. The radial- 
velocity measurements are described in §|3and we determine 
the properties of the system in §01 We compare our results with 
the work by Cocozza et al. (2006) and present the overall dis- 
cussion and conclusions in §|5l In Appendix|X]we elaborate on 
the corrections we applied to the wavelength scale. Finally, we 
discuss the membership of PSR J191 1-5958A with NGC 6752 
in AppendixiBl 

2. Observations and data reduction 

2.1. Spectroscopy 

Twenty-three long-slit spectra of the companion of 
PSR J1911-5958A were obtained with FORSl, the Focal 
Reducer and Low Dispersion Spectrograph of the ESO VLT 
at Cerro Paranal, on 8 different nights from May to August 
of 2004. A summary of the observations is given in Table |2] 
Between the first and second observing run, the instrument was 
moved from Unit Telescope 1 (UTl, Antu) to UT2 (Kueyen). 
The spectra were obtained with the 600 lines mm ' "B" grism 
and a 1'.'31 slit, which gives a wavelength coverage from 
3300-5690A. The standard-resolution collimator was used, 
resulting in a pixel size of 0'.'2pix ' in the spatial direction 
and 1.2Apix ' in the dispersion direction. All spectra had 
integration times of 2470 s and were sandwiched between two 
30 s, B-band, through-the-slit images and preceded by one or 
more 30 s B-band acquisition images. The seeing, as deter- 
mined from the width of the slit profiles, varied between 0'.'5 
and I'.'O, with only three spectra having a FWHM larger than 
C'S. Generally, the conditions were good, with photometric 
skies. With this setup, the spectral resolution is set by the 
seeing, which is less than the slit width in all observations. 
For our average spectra, a resolution of 4.5 A is inferred from 



the spectra of the reference star discussed below. Following 
the FORSl calibration plan, bias, flat-field and wavelength 
calibration frames were obtained during twilight or daytime 
afer each observing night, with the telescope pointed towards 
the zenith. 

Given the proximity of a brighter star {V = 17.3) only 3'.'1 
to the North- West of the pulsar companion, we chose to cen- 
ter the slit on both this star and the companion, see Fig.^ We 
did this to use the star as a reference for the wavelength and 
flux calibration and to minimalize the influence of this star on 
the spectrum of the white dwarf. Besides the pulsar compan- 
ion and this bright star (which we henceforth call the reference 
star or star R), stars A, B and C (see Fig. Q and D also fall 
on the slit. As a result of this setup, the position angle of the 
slit is fixed on the sky and differs from the parallactic angle, 
by an amount which depends on the hour angle of the observa- 
tion. Any effects of differential atmospheric refraction, which 
become important when one does not observe with the slit par- 
allel to the parallactic angle (Filippenko, 1982), are largely cor- 
rected for by the Atmospheric Dispersion Corrector (ADC) on 
the FORS 1 instrument. 

To account for slit losses and to allow for flux calibration, 
the exposure with the slit positioned closest to the parallac- 
tic angle (that from MJD 53229, see Table |2li was followed 
by exposures through a 2" slit of both the pulsar companion 
(160 s) and the spectro- pho tome trie flux standard LTT7987 
(30s: iHamuv et alJl992lll994) . For these the conditions were 
photometric with 0'.'7 seeing. 

The images were reduced with the Munich Image Data 
Analysis System (MIDAS). All images were bias-corrected 
with the bias values from the overscan regions on the FORSl 
chip and flat-fielded using lamp exposures. For the sky sub- 
straction we used clean regions between the stars along the slit. 
The region extended to 26" to each side of the pulsar com- 
panion, encompassing the pulsar companion, the reference star 
and star A, B and C. For star D a similar procedure was used. 
A polynomial was fitted to the spatial profile of these clean re- 
gions of the sky for each column in the disperion direction. The 
order of the polynomial was predominantly zero; but first and 
second order fits were used when this significantly increased 
the goodness of the fit. 

Spectra were extracted from the sky-subtracted image s us- 
ing an optimal extraction method similar to that of H ornd 
(1986). Each of the extracted spectra was wavelength cali- 
brated with the HeHgCd wavelength calibration frames. Here 
we measured the positions of the lines in a row-averaged (in 
the spatial direction) multiplication of the wavelength calibra- 
tion frame and the 2-dimensional (in the spatial and dispersion 
direction) slit profile of the star in question. A cubic polyno- 
mial fit was sufficient to describe the dispersion relation and 
gave root-mean-square (rms) residuals of less than 0.06A. The 
wavelength calibrations were found to be stable between dif- 
ferent observing nights; the rms scatter in the zero point was 
0.041 A (corresponding to 2.7 km s"' at 4500A). 

The last step was to calibrate the spectra for the instru- 
mental response of the spectrograph, as derived from the ob- 
servation of the flux standard. The spectrum of the standard 
was reduced in a similar manner as the pulsar companion and 



C. G. Bassa et al.: The masses of PSR J191 1-5958A and its white dwarf companion 




19 11 46' 19"11'"44=' 19"11"'42* 

Right Ascension (J2000) 

Fig. 1. The orientation of the slit on the sky. The 1'.'31 slit is 
shown with the white-dwarf companion denoted as "WD". The 
nearby star used as a reference is denoted as "R", and names of 
some of the other stars on the slit are indicated. Star D is located 
outside the figure, to the North- West of PSR J1911-5958A. 
This 60" X 60" image is an average of two 30 s B-band acqui- 
sition images obtained during the first observing night. 



the other stars on the slit. Unfortunately, deriving the response 
was somewhat troublesome since the calibrated spectrum of the 
DA white dwarf LTT7987 was tabulated at 50A steps. With 
such a resolution the higher Balmer lines are poorly sampled 
which may result in systematic trends in the flux calibration 
at these wavelengths. We therefore analyzed two archival ob- 
servations of Feige 110 (Qke. 1990), tabulated at lA and 2A 
steps, from June 28th and December 1st, 2004. The observa- 
tions were taken with the same grism as the pulsar compan- 
ion, though with 5" MOS slits. The spectra were extracted and 
wavelength calibrated as before and corrected for atmospheric 
extinction using the average La Silla extinction curve (this re- 
lation is also suitable for Paranal). A comparison of the re- 
sulting response curves showed that these had a very similar 
shape, and that the ratio of the two curves could be well ap- 
proximated with a linear polynomial, i.e. that the response was 
stable over time. We now used the response curve derived from 
the June 2004 observation of Feige 110 and fitted it against 
the extinction-corrected response of LTT 7987, fitting for a lin- 
ear polynomia l scaling factor. Using t he B-band filter curv e of 
iBesselll ( Il990 >) and the zero point of iBessell et"aD h99^ we 
obtain a synthetic B-band magnitude of 12.30 for UrT 7987, 
which compares well to ^=12.2 8 found bv lHamuv etaP jl992l) 
and B=12.27 bv' Landoltllll992i) . 

We corrected all spectra for atmospheric extinction and cal- 
ibrated them using this new response. Differences in continuum 




Fig. 2. Colour-magnitude diagrams of NGC 6752, constructed 
from archival FORSl observations. The stars located on the slit 
are labelled. Shown as solid lines to the left of the cluster main- 
sequence are three Z = 0.0010 helium-core white dwarf cool- 
ing models of Serenelli et al. (2002). The masses of these mod- 
els are as shown (in Mq), and temperatures are indicated along 
the track in panel a, while cooling ages are shown in panel b. 
Also shown are two isochrones from Girardi et al. (2000) for 
an age of 14. 1 Gyr and a metallicity of Z = 0.0004 (solid line) 
and one for an age of 12.6 Gyr with Z - 0.0010 (dashed line). 
All models are placed at a distance of (m - M)v - 13.24 with 
a redd ening of Eb-v = 0.040, as determined bv iGratton et alJ 
(l2003l) . 



flux between the narrow and wide slit exposures were corrected 
for with a scaling factor that depends linearly on wavelength. 

2.2. Photometry 

We have analyzed all available FORS 1 observations of the field 
containing PSR J1911-5958A. The data consists of i) three 
1500 s U, five 360 s B and eight 220 s V-band images, taken 
with the high-resolution collimator (which has O'.'lpix ') on 
3 different nights in 2003 March, April and May under pho- 
tometric conditions with good seeing (0'.'5-0'.'7); ii) three 32 s 
B-band and three 13 s V-band images that were obtained un- 
der photometric conditions with 0'.'7 seeing on 2003 March 3 1 
with the standard collimator; iii) a series of thirty 30 s B-band 
acquisition images obtained prior to the spectral observations 
presented above, and iv) a series of twenty-three 600 s B-band 
images obtained on 2004 August 10-15 with the high reso- 
lution collimator during good to moderate seeing conditions 
(0'.'5-l'.'3). 

All images were bias-subtracted a nd flatfie lded using twi- 
light flats. The DAOPHOT II package llStetsonlll987.) . running 
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Fig. 3. a The B-band magnitude of the white-dwarf compan- 
ion of PSR J1911-5958A as a function of orbital phase for the 
magnitudes determined from the thirty 30 s exposures {open 
squares, offset by -0.2 mag) and the twenty-three 600 s B-band 
exposures {black dots). Also shown are the B-band magnitudes 
determined by ICocozza et alJ ( l2006l) from the same 600 s ex- 
posures {open circles, read off from their Fig. 3, but offset by 
-1-0.7 mag), b The seeing of the 600 s B-band exposures, deter- 
mined from the width of the point-spread-function. Note that 
the variations in the seeing and the magnitudes as measured 
byfCocozza et al. (2006) are highly correlated, indicating that 
their meaurements are coiTupted (see §|5]for details). 



inside MIDAS, was used for the photometry on the a veraged 
images. We followed the recommendations of StetsonI ( fl987.) . 
obtaining instrumental magnitudes through point-spread func- 
tion (PSF) fitting. The B and V-band observations of March 31, 
2003 include 20 photometric standards from St etsoa C20 00) of 
which 12 were unsaturated. The instrumental PSF magnitudes 
of these stars were directly compared against the calibrated val- 
ues to derive zero points and colour terms (no extinction coeffi- 
cients are needed since the standards and target are at the same 
airmass), giving rms residuals of 0.01 mag in B and 0.03 mag 
in V. For the calibration of the U -band observations, we used 4 
standard stars from the PG 1657-H078 field (Landolt, 1992). We 
fitted for zero point and colour term, using the standard ESO 
extinction coefficients of 0.46 and 0.25 mag per airmass for U 
and B-band respectively. These calibrations have rms residuals 
of 0.02 mag in U and B. The difference between the two B- 
band calibrations was less than 0.01 mag, hence we expect our 
zero-point uncertainties in B to be less than 0.02 mag, less than 
0.03 mag in V, and allowing for the uncertainty in the [/-band 
extinction coefficient, less than 0.05 mag in U. The resulting 
magnitudes of the stars on the slit are tabulated in Table [0 



Table 1. VLT/FORSl photometry of the white-dwarf compan- 
ion of PSR J191 1-5958A (denoted with WD) and stars on the 
slit. The nomenclature of the stars is according to Fig.[n while 
star D is located outside the figure, on the North- West side of 
the slit. The uncertainties listed in parentheses are instrumental, 
i.e., they do not include the zero-point uncertainty in the photo- 
metric calibration (0.05 mag in U, 0.02 mag B and 0.03 mag in 
V). The ce lestial positions wer e obtained using the procedure 
outlined in lBassa et alJ ( l2003bl ). 



ID 


Q'2000 


•52000 


U 


B 


V 




h m s 











WD 19 11 42.753 -59 58 26.89 22.02(5) 22.22(3) 22.13(2) 

R 19 11 42.432 -59 58 24.90 17.60(5) 17.78(1) 17.34(1) 

A 19 11 44.768 -59 58 39.58 19.16(5) 19.34(1) 18.80(1) 

B 19 11 42.854 -59 58 33.34 20.36(5) 20.31(1) 19.63(1) 

C 19 11 41.742 -59 58 20.25 18.88(5) 18.84(1) 18.24(1) 

D 19 11 35.911 -59 57 45.35 20.47(5) 20.38(1) 19.70(1) 



Colour-magnitude diagrams were constructed from the 
photometry and are shown in Fig. |2] We find that star R has 
the magnitude and colours of a cluster turn-off star while stars 
A, B and D are located further down the cluster main sequence. 
Star C, on the other hand, is not located on the main sequence, 
and is about a magnitude brighter than cluster stars with the 
same U - V and B - V colours. The pulsar companion is blue 
with respect to the cluster main-sequence by about 2 mag in 
U -V and more than 1 mag in B - V. 

We checked for variability of the pulsar companion using 
the twenty-three 600 s B-band images and the thirty 30 s B- 
band acquisition images. With a pixel scale of CK.'lpix"', the 
600 s images are severely oversampled, and we rebinned the 
images to a pixel scale of 0'.'2pix"' (averaging every 2x2 pix- 
els). Next, instrumental magnitudes were determined through 
PSF fitting and calibrated to the photometry presented above. 
The resulting magnitudes are shown in Fig.|3l In both the 600 s 
and the 30 s images, the magnitudes of the pulsar companion 
do not significantly vary with orbital phase. For example, the 
rms scatter for the white dwarf around the average value in the 
600 s and 30 s images is only 0.02 mag and 0.05 mag, respec- 
tively, and these values are comparable to that of stars of similar 
brightness. 



3. Radial velocities 

Radial velocities of the companion and the five stars on the 
slit were determined by comparing them with template spectra. 
In case of the white dwarf this was done iteratively, where we 
used the hydrogen atmosphere models from which we deter- 
mined the surface gravity and effective temperature (see § 14.21 
and Fig. 13 as a template. A best-fit model was first determined 
for one of the single spectra. This model was then used as a 
velocity template to measure the velocities of the other spec- 
tra. These spectra were shifted to zero velocity and averaged. 
The final velocity template was found by fitting a new atmo- 
sphere model against the averaged spectrum. The actual veloc- 
ities were measured by minimizing a merit function, as de- 
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Table 2. Radial-velocity measurements of the white-dwarf companion of PSR J1911-5958A and four stars on the slit. To put 
these velocities on an absolute scale, a velocity offset of -39 ± 3 km s"' should be added (see AppendixIXt. 











vr 




va 




V'C 




I'D 


MJDb^" 




(km s-') 


(kms-i) 


(km s-') 


(km s-') 


(km s-') 


53147.38143 


0.1590 


-29 ± 22 


12.9 


± 8.1 


18 


± 11 


34 


± 16 


21 


± 20 


53198.11569 


0.7664 


3 ±29 


6.8 


±7.8 


-1 


± 11 


18 


± 16 


11 


± 22 


53204.15036 


0.9752 


-233 ± 26 


11.5 


± 8.2 


23 


± 11 


19 


± 16 


-1 


±20 


53204.19242 


0.0254 


-219 ± 22 


17.8 


± 8.1 


13 


± 12 


45 


± 16 


13 


± 21 


53204.22978 


0.0701 


-160 ± 25 


8.8 


± 8.5 


17 


± 11 


58 


± 15 


12 


± 21 


53206.21868 


0.4459 


218 ±23 


16.8 


± 8.6 


7 


± 11 


33 


± 15 


6 


± 22 


53206.25549 


0.4899 


262 ± 20 


10.8 


± 8.1 


12 


± 11 


38 


± 15 


23 


± 21 


53206.29500 


0.5371 


245 ± 23 


10.3 


± 8.5 


3 


± 11 


31 


± 16 


10 


± 20 


53210.19705 


0.1983 


6 ±25 


15.4 


± 8.8 


2 


± 11 


35 


± 16 


20 


± 20 


53210.23098 


0.2388 


12 ±22 


-0.3 


± 8.0 


-1 


± 12 


30 


± 15 


3 


± 20 


53210.27816 


0.2952 


77 ±24 


-6.4 


±7.6 


12 


± 11 


38 


± 16 


2 


± 21 


53210.31250 


0.3362 


148 ± 27 


3.6 


±7.8 


-1 


± 11 


33 


± 16 


23 


±26 


53229.01322 


0.6748 


116 ±23 


6.2 


±8.1 


8 


± 11 


27 


± 16 


-1 


± 22 


53231.02712 


0.0805 


-186 ±25 


0.3 


±8.1 


8 


± 11 


12 


± 16 


-12 


± 22 


53231.06463 


0.1253 


-179 ±27 


3.8 


± 8.2 


6 


± 11 


19 


± 16 


-8 


± 20 


53231.11869 


0.1899 


-64 ± 26 


16.3 


±7.8 


6 


± 11 


33 


± 16 


22 


± 21 


53231.17422 


0.2562 


35 ±25 


-7.0 


±7.9 


4 


± 11 


22 


± 16 


-6 


±20 


53231.22574 


0.3177 


153 ± 29 


7.7 


± 8.2 


-6 


± 10 


34 


± 16 


-6 


± 21 


53232.00988 


0.2544 


57 ± 22 


1.2 


±7.7 


5 


± 11 


9 


± 16 


-8 


± 19 


53232.04617 


0.2978 


42 ±24 


11.4 


±8.1 


13 


± 11 


21 


± 15 


9 


± 20 


53232.08482 


0.3439 


142 ± 21 


23.9 


± 8.2 


15 


± 10 


38 


± 15 


16 


± 20 


53232.13559 


0.4046 


200 ± 22 


12.7 


± 8.4 





± 11 


37 


± 17 


-5 


±20 


53232.18610 


0.4649 


273 ± 32 


-0.6 


±9.2 


3 


± 10 


32 


± 17 


-1 


± 21 



" The time of the observ ation at mid-exposure, corrected to the solar system barycenter. 

* Using the ephemeris of iD' Amico etall <2002h : T.^c = MJD 51919.2064780(3), Pt = 0.837113476(1) days. 



fined in lBassa et al ] (l2003ah . fitting for velocity and a 2nd order 
polynomial modelling continuum differences. 

In the case of the reference star we used a somewhat dif- 
ferent approach. Here, a template was constructed from fitting 
Lorentzian profiles to a single, normalized spectrum of the ref- 
erence star. Eight lines (H/3 up to Hll and Ca K, but without 
the blend of He and Ca H) were simultaneously fitted, fitting 
for depth and width, but forcing the velocity to be the same for 
all lines. The resulting template was shifted to zero velocity and 
fitted against normalized spectra of the reference star by again 
minimizing a merit function. To test for the stability of the 
radial velocities, we also used this template to determine radial 
velocities of the other four stars on the slit. These stars also dis- 
play the hydrogen Balmer lines and Ca K, though the lines are 
not as strong as those of the reference star (see Fig.|^3. 

The radial velocity of the pulsar companion varied by 
as much as 470km s"' between the different spectra, fully 
within the expectations for this system. From radio timing, 
it was found th at the eccentricity of the orbit is e < 10"^ 
jP' Amico et al.l |^02). Hence, we fitted the barycentric ra- 
dial velocities of the pulsar companion to a circular orbit, with 
the orbital period and time of the ascending passage node 
fixed t o the values de t ermine d from the radio-timing solu- 
tion by ' D' Amico et al] (l2002l) . We find a radial-velocity am- 
plitude /TwD - 215 ± 10km s~' and a systemic velocity y - 
-30 + 7kms"' for a reduced chi-squared = 1.67 with 21 



degrees-of-freedom. Here, the errors on A^wd and y are scaled 
to give = 1. 

The velocity of the reference star, however, varied over a 
120 kms"' range, much larger than the 7 to lOkms ' error on 
the individual velocities. We feared that these velocity varia- 
tions might also be a result of binarity, but the other stars on 
the slit displayed similar variations in velocity. In particular, 
the velocities of these stars displayed a trend when compared 
against the local hour angle of the observation, where the ve- 
locity decreased by about 16 km s ' per hour prior or after cul- 
mination. Upon closer investigation this trend was found to be 
caused by two separate, systematic, effects. Because the effects 
are systematic, they can be corrected for. 

For the first correction, we applied a shift to the wavelength 
calibration of each spectrum based on the difference between 
the measured wavelength of the O i A5511 night sky emission 
line in the spectrum and the labaratory value. This correction 
removes the global decrease of the velocities as a function of 
hour angle. The remaining scatter in the radial velocities is 
largely removed by correcting for the second effect, which is 
due to eiTors in the centering of the stars on the slit. For this 
correction, we determined the position of star R with respect 
to the center of the slit in the through-the-slit images obtained 
before and after each spectrum. We apply this offset as a shift 
in wavelength to the wavelength calibration of the correspond- 
ing spectrum. In AppendixlAlwe describe these corrections in 
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Orbital Phase 

Fig. 4. The radial velocities of the white-dwarf companion of 
PSR J1911-5958A {black dots) and of the reference star R 
(open circles). The solid line represents the best-fit model for 
a circular orbit, using all data points, with the time of the as- 
cending node passage and the orbital period fix ed to the values 
determ ined from the radio-timing solution by D'Amico et alJ 
ll2002 h. The dashed line represents the best-fit model exclud- 
ing the boxed data point. To put these velocities on an absolute 
scale, a velocity off'set of -39 + 3 km s"' should be added (see 
AppendixIXt. 

detail, while Tabled lists the velocities that were determined 
from the corrected wavelength calibrations. 

We now use the corrected radial velocities to determine 
the radial-velocity orbit of the white-dwarf companion of 
PSR J1911-5958A. Again fitting for a circular orbit, we find 
a radial-velocity amplitude /^wd - 226 + 9kms"' and a sys- 
temic velocity y = 26 + 6kms"' (xl - 1-45 for 21 degrees- 
of-freedom). The errors on /Two and y are again scaled to give 
xl = I- This fit is represented in Fig. 0] with the solid line. If 
we exclude the single point (the boxed point in Fig.|4} that lies 
3.2cr away from the best-fit, the fit improves to xl - 1 -00, giv- 
ing /TwD - 231 ± 8kms"' and y = 21 + 5kms"', depicted 
by the dashed curve in Fig. |31 Remarkably, this outlier corre- 
sponds to the spectrum taken during the first observing run, 
when, as mentioned in §|5] FORSl was still on UTl, unlike 
all other spectroscopic observations, when it was on UT2. We 
do not understand, however, how this could cause a difference, 
since differences in the flat-fields, wavelength calibration, and 
flux calibration should all be corrected for. 

If we fit a circular orbit against the radial velocities of the 
pulsar companion relative to the radial velocities of the refer- 
ence star R, we obtain /Tav = 225 + lOkms, 7ai' - 17±6kms"', 
xl = 1.43 for 21 degrees-of-freedom. Again excluding the out- 
lier gives K/^y = 231 + 9kms"\ -/a,, = 13 ± 6kms"', with 
xl = 1.10 for 20 degrees-of-freedom. 



The differences in the radial velocity amplitudes /Two are 
consistent within the errors. The same holds for the systemic 
velocities y measured from the absolute and the relative ve- 
locities. For the remainder of this paper we will use the fit 
using absolute velocities, without the 3.2cr outlier; /Two - 
231 ± Skms"' and y - -18 ± 6kms^' (here we corrected the 
systematic velocity for the -39 ± 3 km s"' velocity off'set which 
we determined in AppendixlAl. 



4. System properties 

We use our measurements to determine the properties of the 
white dwarf and the pulsar In our analysis, we distinguish be- 
tween results that do and that do not depend on the assumption 
that the pulsar is a member of NGC 6752. We will see that our 
conclusions depend on that assumption; we will address this in 
detail in AppendixiBl 



4.1. Minimum wiiite dwarf mass 



The radio timing observations by b'Amico et all (l2002h 
yielded a projected semi-major axis of the pulsar orbit of 
flpsR sin i/c = 1 .206045 + 0.000002 s, which, together with 
the orbital period Pb implies a radial-velocity amplitude of 
KpsR = 31.40986 ± 0.00005 km s"'. Combining this with the 
radial-velocity amplitude of the white dwarf determines the 
mass ratio q = Mpsr/Mwd = -KVo/^psr = 7.36 + 0.25. Here, 
the error on q is dominated by the uncertainty in /Two- 

We can use the mass ratio and the constraint that the incli- 
nation must be less than or equal to 90° to determine a lower 
limit to the white-dwarf mass. For this, we use the pulsar mass- 
function /(Mpsr) = sin^ //(Mwd + Mpsr)^ = (2.687603 + 
0.000013) X 10"^ Mq, so that we can write Mwd sin"* / = (1 + 
^)^/(Mpsr). Setting the inclination at ; - 90° and using the 
mass ratio q as determined above we find a Icr lower limit of 
Mwd > 0.177 Mq. The 2cr lower Hmit is Mwd > 0.166 Mq. 

4.2. Effective temperature and surface gravity 

The atmospheric parameters for the white dwarf were deter- 
mined by fitting theoretical model atmospheres to the aver- 
age of the velocity-corrected spectra (see § O. The theoret- 
ical models were taken from a grid of pure hydrogen mod- 
els usually applied to normal DA white dwarfs, but extending 
down to surface gravities of logg - 5. The methods and i nput 
physics are des cribed in more detail in ' Finlev etalJ ( ll997l) and 
iHomeier et alJ 0199 8). The best-fitting parameter s are found 
with a Levenberg-Marquardt type x^ algorithm jPress et all 
'1992^. We also use a second, newly developed algorithm, 
which is less sophisticated but more transparent and robust than 
the Levenberg-Marquardt method. In essence it determines the 
X^ values for the models of the grid around the minimum and 
then fits the x^ surface with a paraboloid, from which the pa- 
rameters and errors corresponding to the minimal can be cal- 
culated. This avoids excessive interpolations between the mod- 
els in the grid, which sometimes leads to artificial small-scale 
structure of the;if^ surface. The results between the two meth- 
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ods did not differ significantly; the values given below are from 
the second method. 

For the fit we used the spectral ranges from 3740^440 and 
4760-5030 A, which contain the Balmer lines. The model was 
fitted to the observed spectrum and the calculated from the 
fit to the continuum-normalized line profiles. The obtained pa- 
rameters and their formal errors are Teff = 10090 + 25K and 
\ogg = 6.44 + 0.05 cgs. The resulting best-fit model is shown 
in Fig. |5l With this temperature, the spectral type of the white 
dwarf is DA5 ( Wesemael et aT . 1993). 

The model fits the observed Balmer lines extremely well 
up to HIO. Hll and H12 are weaker and not as deep in the 
model as is observed. This may be an indication that the sur- 
face gravity is slightly lower than the formal fit. It may also in- 
dicate that the Hum mer-Mihalas occupatio n probability theory 
llHummer & Mihalas..l988:.Mihalas et al.Lll988„ 1990J used in 
the models overestimates the quenching of the highest Balmer 
lines. The theoretical line profiles of Hll and H12 depend 
strongly on the Hummer-Mihalas formalism and as these lines 
are absent in the typical logg - 8 white dwarf, the theory 
is difficult to calibrate. For a recent study on the influence 
of the Hummer-Mi halas formalism on the line profiles, see 
iKoesteretaLllEooi. 

The fit uncertainties include only the stati stical errors, 
while system atic ef fects may be important feeidl Il996t 
[Vauclair et al. I Il997l) . We experimented with small changes 
in the continuum and fit regions, and other fitting parame- 
ters. The largest effect is from changes in the resolution: as- 
suming that the resolution were 6.5A instead of the 4.5A in- 
ferred from the lines in the spectrum of the reference star, we 
find Tetf = 10135 K and logg = 6.20cgs. As we discuss in 
Appendix IB. 2. 21 rapid rotation may mimic a change in reso- 
lution. Since this is not secure, however, we decided to adopt 
conservative errors below, and use Teff = 10090 ± 150K and 
logg = 6.44 + 0.20 cgs 

4.3. White dwarf mass and radius 

The mass and radius of the white dwarf can be determined 
from the surface gravity using a mass-radius relation appro- 
priate for a helium-core white dwarf at the observed tem- 
perature. Such m ass-radius relations have been modelled by 
IPanei et a r(2000) and we interpolate between their 8000 K and 
12000K tracks to obtain the relation at the observed temper- 
ature, as shown in Fig. |6l giving Mwd = 0.180 ± 0.018 Mq 
and 7? WD = 0.0423 ± 0.0075 Rq. Similar white-dwarf masses 
and radii are found from the evolutionary cooling tracks of 
[Oriebe et al. (1998.1 and iRohrmann et al. (2002). Here we ob- 
tain, for each model with a given mass, the radius and hence the 
surface gravity at the observed white-dwarf temperature and 
interpolate between the models t o get the mass and radi us at 
the observed log g. The models bv lRohrma nn et ar(2002) give 
Mwd = 0. 1 8 1 ± 0.0 1 2 Mo with Rwr , ^ 0.0 424 ±0.0088 Rr,. The 
lowest mass model bv lDriebe et al. I ll 19981) has 0. 1 79 Mq, so we 
extrapolate their models. This yields Mwd = 0.172 + 0.018 Mq 
and /?WD = 0.0414 ± 0.0074 7?o, though the uncertainties may 
be underestimated because of the extrapolation. The uncer- 



tainties of the lRohrmann et aD (l2002l) models are considerably 
smaller, as i ts mas s-ra dius relation is steep er than those of the 
'Panei et al.' ('2000') and lPriebe elaP (Il998l) . Finafly, the mass- 
radius relations from the Z = 0.0010 and Z - 0.0002 models by 
ISerenelhet alJ l l2o61 are very similar and both predict a some- 
what higher white-dwarf mass, of Mwd = 0.190 + 0.015 Mq, 

and a radius of /?wd = 0.0434 ± O.OO84£0; 

We should note tha t the models bv lDriebe et all ll 19981) and 
iRohrmann et al ](|2002') are computed for white-dwarf progen- 
itors with solar metallicity and are appropriate for field sys- 
tems. As the metalhcit y of NGC 6752 is c onsi derably smaller 
([Fe/H]= -1.43 + 0.04, Gratton et al.'2003'), the' SereneUi etail 
(2002) models for white-dwarf progenitors with sub-solar 
metallicities (with Z = 0.0010 in particular) are more appro- 
priate in the case that PSR J1911-5958A is associated with 
NGC 6752. The diff'erences of about 0.01 Mq and 0.001 Rq in 
the predictions from different models, however, are similar in 
magnitude to the difference induced by different metallicities. 
For the remainder of the paper, we will use values that encom- 
pass all predictions from the effective temperature and surface 
gravity: Mwd = 0.18 + 0.02 Mq and /?wd = 0.043 + 0.009 Rq. 

4.4. Distance inferred from the white dwarf 

The distance to the white dwarf can now be estimated using the 
observed and modelled flux and the radius of the white dwarf. 
We do this by writing the flux normalization / = {R/d)^nF 
between the observed flux /, the model flux nF, and the ra- 
dius R over the distance d in terms of magnitudes; M^ - 
43.234 - 5\ogR/RQ - l.SlognF^ + c,i, where M^ is the ab- 
solute magnitude in a given filter, ttF^ the flux from the model 
in the same filter and the zero-point of the filter. By con- 
volving the flux-calibrated best-fit model of the observed spec- 
trum with the B and V-band filter curves from .BesselL il990l) . 
we obtain ttFb - 9.45 x 10^ ergcm"^ s"' A""' and nFy = 
6.26 X 10^ ergcm"^ s"' A"'. The uncertainties on these fluxes 
due to the uncertainties in Tetf and logg are about 5%. With 
the zero-points from Bessell et al. ( 1998), cb = -20.948 and 
cv = -21.100, we obtain Mb = 9.63 ± 0.46, My = 9.48 ± 0.46. 
The uncertainties are dominated by those on the white-dwarf 
radius. 

Combined with the observed B and V-band magnitudes 
from Table we obtain distance moduli of (m - M)b - 
12.59 + 0.46 and (m-M)v = 12.66 + 0.46. Assuming a negligi- 
ble reddening, the averaged distance modulus yields a distance 
of li = 3.4 ± 0.7 kpc. Assuming a reddening of Eb-v = 0.05, 
the distance becomes d = 3.1 + 0.7 kpc. 

4.5. Mass and radius from distance, flux and 
temperature 

Now we determine the properties of the white dwarf based on 
the assumption that the distance of the binary is that of the glob- 
ular cluster. From the values from Tabled the distance modu- 
lus (m - M)v ^ 13. 24 ± 0.08 and reddening Eb-v = 0.046 ± 
0.005 jGratton et all. 2003. ■2005.) . we obtain Mb = 8.93 ± 0.09 
and Mv = 8.89 + 0.08. The B and V-band fluxes and zero- 
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Fig.5. Spectrum of the white-dwarf companion to PSR J1911-5958A. The lower curve in panel a shows the average of the 
23 individual spectra, shifted to zero velocity. The hydrogen Balmer lines are seen from H/3 up to H12 as indicated. The top 
curve is the best-fit model spectrum, with T^ff - 10090K and logg - 6.44 cgs. It is shifted upwards by 5/iJy. Panel b shows the 
flux-normalized line profiles superposed on the modelled profiles. The Hfi profile is shifted a quarter unit downwards from unity, 
while those of H6 up to H12 are shifted upwards by multiples of the same amount. 



Driebe et al. 1998 
Panei et al. 2000 
Rohrmann et cl. 2002 
Serenelli et al. 2002 




Fig. 6. White-dwarf mass-radius relations for a temperature 
of Teff - 100 9 + 1 50 K. Shown a re the relation s frorn 
iRohrmann etaP (l2002i) (dashed lin e), ISereneUi et aP (l2002h 
(dashe d dotted), IPriebe et"^ lll998l) (dotted) and IPanei et alJ 
hoOO*) (solid lines). For the latter model, the 8 000K and 
12000K models used to obtain the mass-radius relation at the 
observed temperature are also shown. The diagonal solid and 
dashed lines depict the observed range (Icr) in surface gravity 
(logg = 6.44 + 0.20 cgs, with g = GM/R^). The horizontal 
solid and dashed lines indicate the white-dwarf radius deter- 
mined using the distance of NGC 6752. The filled light grey 
area depicts the region excluded by the 2cr lower limit on the 
white-dwarf mass (Mwd > 0.166 M©) that was derived from 
the pulsar mass-function and the observed mass ratio. 



points and the relation from § 14.41 yield white-dwarf radii of 
\ogR/RQ = -1.226±0.021 and logR/R^ = -1.249 ±0.018 for 
the B and V-band values, respectively. As a conservative esti- 
mate, we will use log R/Rq = -1.238 + 0.030, corresponding 
to R = 0.058 + 0.004 7?Q, which encompasses both values. 

The mass of the white dwarf can be determined from the ra- 
dius using the mass-radius relations. The radius is shown with 
the horizontal lines in Fig.|6l together with the mass-radius re- 
lations. The models bv Rohrmann et al.. (2002) predict Mwn ^ 
0.172 ± 0.001 Mo, while the models of ISereneUi eTaP ll2002l) 
give Mwn = 0.175 ± 0.002 Mq. The mass-radius relatio n by 
IPanei et al.. (2000) and the evolutionary models bv Driebe et alJ 



(11998J) do not reach these radii and provide no mass estimate. 
The surface gravity that coiTesponds with these masses and 
radii is about log g = 6.20 cgs. Both the white-dwarf radius and 
the surface gravity are slightly outside the Icr range of these 
values inferred from the spectrum (see Fig.|6j. Because of the 
steepness of the mass-radius relations, the masses are in agree- 
ment. We do note that if the pulsar binary is at the distance of 
NGC 6752, the white-dwarf mass is slightly below our best es- 
timate for the minimum mass, though they are consistent with 
the 2cr limit. 

If the pulsar binary is associated with NGC 6752, the pro- 
genitor of the white dwarf m ust have evolved from a clus- 
ter star. ISereneUi et alJ (l2002l) computed white-dwarf cooling 
tracks for objects with sub-solar metallicities. Their Z = 0.001 
models have a metallicit y that is similar to that of NGC 6752 
([Fe/H]= -1.43 ± 0.04; Gratton et alJl2003l) . and for this rea- 
son, we use the coiTesponding mass estimate as the best value 
for the white-dwarf mass. To take into account the uncertainty 
in the white-dwarf mass-radius relation, we add 0.010 in 
quadrature to the uncertainty in the mass. In summary, un- 
der the assumption that the binary system is a member of 
NGC 6752, we infer a radius of /?wd = 0.058 ± 0.004 /^q and a 
mass of Mwd = 0.175 + 0.010 Mq. 



C. G. Bassa et al.: The masses of PSR J191 1-5958A and its white dwarf companion 



9 




Fig. 7. The various constraints on the mass of 
PSR J1911-5958A. Panel a shows the constraint set by 
the mass ratio as the diagonal solid and dashed black lines. The 
large filled light grey area is excluded by the mass-function 
of the pulsar, as points in this area would require sin / > 1 . 
The line that borders this area is for the limit / - 90° in the 
pulsar mass-function. The two dotted lines depict inclinations 
of / = 75° (lower line) and / - 60° (upper line). Two pairs of 
vertical error bars on the right-hand side of the panel represent 
the two white-dwarf mass estimates and their uncertainties, 
where the larger error bar is for the mass determined from the 
surface gravity and eff'ective temperature, whereas the smaller 
error bar is from the assumption that the binary is a member of 
NGC 6752. Allowed values for the pulsar mass exist in the light 
grey region for a white dwarf associated with NGC 6752, and 
the light and dark grey area for a white dwarf not associated 
with the globular cluster The resulting pulsar masses, based on 
the Monte Carlo simulation, are indicated with the horizontal 
error bars. The uncertainties on the white-dwarf mass, pulsar 
mass, mass ratio and the mass-function are all 95% confidence 
(2cr). The probability densities as a function of pulsar mass are 
shown in panel b. The 68% and 95% confidence uncertainties 
based on these distributions are denoted by the open squares 
and the filled circles, respectively. 

4.6. Pulsar mass 

As for the white dwarf, we can determine a minimum mass 
for the pulsar from the mass ratio and the constraint that the 
inclination / is less than or equal to 90°. The white-dwarf 
mass-function is /(Mwd) = ■'^psr ^i^^^ (-^wd + A/psr)^ = 



^wD^b/(27rG) = 1.072 + 0.108 Mq, which can be written to 
give MpsR sin-' / = (1 + 1 /^)V(^wd)- With sin / < 1, this yields 
a Icr lower limit to the pulsar mass: MpsR > 1.24 Mq. The 2cr 
limit is MpsR > 1.10 Mq. These limits are model independent, 
as they only depend on three observables; the projected semi- 
major axis of the pulsar orbit apsp sin / and the orbital P^, as 
determined from pulsar timing, and the white-dwarf velocity 
amplitude /Two, that we obtained from our spectroscopic ob- 
servations. 

In Fig.0 we show the constraints set by the mass ratio and 
the two white-dwarf mass determinations. It is clear that a large 
part of the range allowed by these constraints is excluded by 
the pulsar mass function. As a result, the most probable value 
for the pulsar mass and the uncertainties on it are not normally 
distributed. Instead, we determined these values via a Monte 
Carlo error propagation method. For a million trial evaluations, 
values for P^, apsR sin /, /Two and Mwd were randomly drawn 
from Gaussian distributions with the appropriate means and 
widths to obtain the corresponding pulsar mass. Solutions that 
had sin / > 1 were excluded. From the resulting distribution of 
solutions the most probable value for and the uncertainties on 
the pulsar mass were determined. 

For a white-dwarf mass of Mwd = 0.18 ± 0.02 Mq, the 
mass of the pulsar is constrained to Mpsp - 1.4O^q[qM0 at 
68% confidence. For the case that the white dwarf is associated 
with the globular cluster, the allowed range in pulsar mass is 
smaller, 1.34 + 0.08 Mq at 68% confidence. The uncertainties 
corresponding to 95% confidence are in both cases twice as 
large. 

5. Discussion and conclusions 

We have unambiguously identified the companion to 
PSR J1911-5958A as a Helium-core white dwarf and 
determined its mass. Together with the measurement of the 
mass-ratio of the binary, we obtain constraints on the pulsar 
mass. However, before we d iscuss our results, we compare our 
results to those presented bv' Cocozza et alJ ( l2006h . 

Cocozza et al. found that the light-curve of the white dwarf 
companion to PSR J 191 1-5958 A showed two phases of bright- 
ening by about 0.3 mag, located approximately at the quadra- 
tures of the orbit (phases = 0.0 and = 0.5). This result is 
at odds with our light-curve, which excludes variations larger 
than 0.02 mag. Fig.E] shows a repr oduction of the light-curve 
determined bv ICocozza et al.l (1200 6). This figure also shows 
the variation in the seeing under which these images were ob- 
tained; one sees that these closely follow the var i ation s in the 
white-dwarf magnitude found by Cocozza et al.l (1200 6). This 
suggests that seeing aff'ects their photometry and that the vari- 
ations they measure are due to variations in the seeing and not 
due to variations in the white-dwarf brightness. We found that 
we could reproduce their light-curve by defining the PSF over 
an area smaller than about 3'.'1 in radius. This distance corre- 
sponds to the distance between star R and the white dwarf, and 
if the PSF radius is chosen smaller than this value, flux in the 
wings of star R is not removed and added to the flux of the 
white dwarf. Hence, the effect increases for increasing seeing. 
We used a PSF radius that extends up to 4" from the center of 
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each star and is still 2.5x larger than the width of the PSF in 
the images of the worst seeing. As such, our photometry is not 
affected by this error and excludes the 0.3 mag variations seen 
by Cocozza et al. (2006). 

The radial-velocity curve of the whi te dwarf companion t o 
PSR J1911-5958A is also presented in lCocozza eTal] ilOQ^ . 
Their radial-velocity amplitude /Two and systemic velocity 
■y are consistent with the values we found using the uncor- 
rected velocities. We note that the uncertainty on our value for 
/TwD is about a factor two smaller. We believe this is caused 
by the f act that we used nine Balmer lines (H/3 up to HI 2) 
whereas fcocozza et alJ ilQOd ) only used four (H/3 up to He). 
Especially since the higher Balmer lines are narrower they 
will provide more accurate velocities. lCocozza et al.l (2006) use 
their measurement of the systemic velocity of the pulsar binary 
(y = -28.1 ± 4.9kms"') as an arguement supporting the as- 
sociation of PSR J19 1 1-5958A with NGC6752. However, as 
ICocozza et alJ ll200 (j) did not correct for the systematic shifts 
in the wavelength scale that we identified and corrected for 
(see Appendix O, their conclusion regarding the association 
between the binary and the cluster is meaningless. 

We now turn to the conclusions that can be drawn from our 
results. In AppendixlEl we have used the available constraints 
set by our observations to test whether PSR J1911-5958A is 
associated with the globular cluster NGC6752. Unfortunately, 
these tests are not conclusive and hence, we discuss both pulsar 
mass determinations below. 

First, it is interesting to compare the mass of the white 
dwarf with the mass predicted by the theoretical r elation be- 
tween the white-dwarf mass and the orbital period Jjoss et all 
Il987h . For short orbital periods, this relation is least secure, 
since mass tran sfer sta rts before the companion has a well- 
developed core ( lErgma et aL.199 8). Nevertheless, from earlier 
systems it seemed that the predictions by Tauris & Savoniie 
lll999h . which are strictly valid only for > 2d, work 
well for binaries with orbital per iods as short as 8h (see 
Fig. 2 in Ivan KerkwiiketallbOOSh . At the orbital period of 
PSR J1911-5958A, their relation predicts a white-dwarf mass 
between 0.18 and 0.20 Mq. This is again in very good agree- 
ment with our white-dwarf mass measurement (independent of 
whether the pulsar binary is associated with the globular clus- 
ter). 

Binary evolution furthermore predicts that the progenitors 
of white dwarfs in low-mass binary millisecond pulsars have 
lost ^ 0.6 Mq of matter in order to form a ~ 0.2 Mq helium-core 
white dwarf. It is believed that at least a part of this matter is ac- 
creted onto the neutron star in order to spin it up to millisecond 
periods. As such, the neutron stars in low-mass binary pulsar 
systems are expected t o be heavier than the canonical neutron 
star of 1 .35 ± 0.04 Mq (Tho rsett & Chakrabartvlfl99^ . 

For the case that PSR J1911-5958A is a field sys- 



tem, the mass of the pulsar (MpsR = 1.40+^;|°Mo) IS m- 
deed heavier than the canonical value, though not by much. 
However, similarly small differences have been found for 
several of the other low-mass binary millisecond pulsars 
for which masse s have bee n determined; PSR J1713+0747, 
1.3 + 0.2 t.Splaveret ali>2005 ). PSR J1909-3744, 1.438 ± 
0.024 Mo iTacobv et all l2005l) . PSR J0437-4715 with 1.58 + 



0.18 dvanStraten eta4 1200 ih . PSR B1855H-09, 1.6 + 
0.2 Mm jNiceet al.. 2005a) PSR J1012-H5307, 1 . 6 + .2 Mq 
Jvan Kerkwiik et aU il99a .2001 ICallanan et all Il998l) . The 
only system for which the pulsar is significantly h eavier than 
the 1.3 5 Mo is PSR J075 1 -i- 1 807, with 2. 1 + 0.2 Mq jNice et all 
l2005bl) . 

For the case that PSR J1911-5958A is associated with 
NGC 6752, the pulsar mass ( 1 . 34 + 0.08 Mq) is consiste r it with 
the 1.35 + 0.04 Mo found bv iThorsett & Chakrabartvl ( 1 19991) 
and is one of the least heavy neutron stars in low-mass binary 
millisecond pulsars. In this case, it is interesting to compare 
PSR J1911-5958A with PSR J0737-3039B, the 2.8 s non- 
recycled pulsar in the double pulsar system. If the mass of this 
pulsar (1.250 ± 0.005 M^i lLvne et alJ2004 is indicative of the 
mass of a neutron star after is has been formed, it would only 
take less then 0. 1 Mq to recycle a normal neutron star to a mil- 
lisecond pulsar spinning with a period of 3.26 ms. 

Finally, our observations also constrain the inclination of 
the system. For the case that PSR J1911-5958A is not asso- 
ciated with NGC 6752, we have a 2cr limit of sin / > 0.923 or 
/ > 67?4. In the other case, sin / > 0.968 and / > 75^5. Because 
of these high inclinations, the effects of Shapiro delay should 
be significant in the timing of the pulsar. Combined with our 
measurement of the white-dwarf mass, these limits on the in- 
clination imply a Shapiro delay As > 5.7 fj.s. Unfortunately, for 
nearly circular orbits, the Shapiro delay is covariant with the 
projected semi-major axis and the eccentricity, and a large part 
of the effect is absorbed in these two parameters. As a result, 
the effect that remains has a size of Aj' > 1.2 fis. Interestingly, 
if Shapiro delay is present, but neglected in the pulsar timing 
fit, it introduces an apparent eccentricity of e > 1.3 x 10"^ and 
places the longitude of periastron ato) - 90°. Though the small 
signal due to Shapiro delay may be difficult to detect, radio tim- 
ing observations of PSR J191 1-5958A may be used to obtain 
an upper limit on the inclination and the companion mass. For 
example, if / - 85°, the timing signal due to Shapiro delay will 
be much larger. As - 11.2/is and As' = 5.4 //s. 
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Appendix A: Corrections to the wavelength 
calibration 

First, according to the wavelength calibration of the different 
spectra, the wavelength of the Oi /15577 night sky emission 
line was offset from the laboratory value (A = 5577.34 A). 
The offsets varied over a range of 0.9A or about 50kms"' and 
appeared to decrease with increasing hour angle, as shown in 
Fig. lA.lb . We believe this effect is the result of instrument flex- 
ure due to differences in the pointing of the telescope, as the 
calibration frames are obtained with the telescope pointing to- 
wards the zenit, while the PSR J191 1-5958A has secz > 1.24. 
We corrected for this effect by applying the measured offsets as 
a wavelength shift in the zero point of the wavelength calibra- 
tions of each individual spectrum. 

However, as the O i /15577 sky line lies redwards of the last 
arc line used for the wavelength calibration (Hg /15460) we 
cannot verify that the polynomial used for the wavelength cal- 
ibration is still valid at the wavelength of the skyline and that 
all systematic wavelength shifts are corrected for. We found, 
however, that the higher order terms of the polynomial are very 
similar between the calibration spectra extracted for the differ- 
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ent stars; this is because the curvature of the cahbration Hnes 
and the sky line at the positions of the stars is very small. As a 
result, we do not expect random effects to dominate the wave- 
length shift computed above, though the procedure may intro- 
duce a systematic wavelength shift affecting all spectra. 

The second effect involved the centering of the stars on the 
1'.'31 slit. In general, when a star is not centered on the slit the 
offset of the star from the center of the slit will appear as a 
wavelength shift in the spectrum. This can clearly be seen by 
the velocities of star R and A in Fig. lA.lh . where the velocity 
of star A differs from that of star R by about 100 km s"' ; this is 
primarily due to the fact that compared to star R, star A is po- 
sitioned closer to the right-hand (redwards) edge of the slit, see 
Fig-E Because of the magnitude of this effect, we determined 
the position of the reference star R with respect to the center 
of the slit from the through-the-slit images taken before and 
after the spectra. The positions on the 'before' images varied 
over a range of 0.68 pixels, whereas those on the 'after' im- 
ages varied over 1.81 pixels. Interestingly, the 'after' positions 
with negative hour angles displayed positive shifts (redwards in 
wavelength) whereas those with positive hour angles displayed 
negative shifts. Fig. IA.lt shows these trends. We interpret the 
variations in centering of the reference star R on the 'before' 
images as simple scatter inherent to the centering of a star on 
the slit. The centering variations on the 'after' images clearly 
has a different cause, given its dependence on the hour angle. It 
may be that this is also related to flexure or differential atmo- 
spheric refraction. The averaged affect introduces the scatter in 
the velocities of the stars on the slit (Fig. lA.lh ) and affects all 
stars in the same way. 

To correct for this effect, we computed the position x, of 
star i (in pixels) with respect to the center of the slit and applied 
it as a shift in wavelength in the zeropoint of the wavelength 
calibrations of each individual spectrum. The position is com- 
puted from Xi - |/'psf(-^- -^R - ^■^z) djc, where (/^psf is the 
normalized point-spread function, as determined from fitting a 
Moffat function to the spatial stellar profile. Furthermore, jcr is 
the average of the 'before' and 'after' shifts in the position of 
the reference star R with respect to the center of the slit and w 
is the slit width (in pixels). Finally, Ax, is the offset between the 
position of the reference star R and star /, measured in pixels in 
the dispersion direction. These offsets were computed from the 
five combined 360 s B-band exposures used for the photometry. 

The resulting radial velocities of the pulsar companion and 
the other stars on the slit, corrected to the solar-system bary cen- 
ter, are given in Table |2] The velocity of the reference star R 
now only varies over a range of 30km s ', with an rms scatter 
of 11 kms"', comparable to the errors on the velocities tabu- 
lated in Table |2l For stars A, B, C, and D the rms scatter is 
comparable or somewhat larger, with 10, 50, 32 and 13 kms"', 
respectively. The large scatter in the radial velocities of star B 
is likely caused by secondary light from the two nearby stars 
(see Fig.Q, where variations in the seeing lead to shifts in the 
center of light in the dispersion direction, which in turn leads to 
velocity shifts. Therefore, we have not used star B in the further 
analysis. 



T ' 1 ' 1 ' 1 ' 1 ' r 
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-2-101234 

Hour angle (h) 

Fig. A.l. The systematic effects that were present in the radial 
velocity study, a The uncorrected radial velocities of stars R 
and A as a function of hour angle. For both stars, the velocities 
decrease with increasing hour angle, while the scatter around 
the main trend is comparable for both stars, b The offset in 
the measured position of the O i A5577 sky line compared to 
the laboratory value, c The centering of star R on the slit, as 
measured from the 'before' and 'after' through-the-slit images. 
Each line connecting the two black dots correspond to a single 
spectrum, where the left dot is the position determined in the 
'before' image, and the right-hand dot the position in the 'after' 
image. The open circle is the average of the two positions in 
hour angle and pixel shift. 

The stars on the slit have systemic radial velocities of 
yR = 8 + 2kms"', ja - T ± 2kms"', yc = 30 ± 2 km s"' and 
yo =6 + 2 km s"' . It is unexpected that all these stars have sys- 
temic velocities that are different from the radial velo city of the 
globular cluster NGC6752 (-32.0 ± 1.6kms'': iDubath et alJ 
1997), especially as all stars, except star C, have colours that 
place them on the naiTow cluster main sequence in Fig.|2] Since 
star R, A and D have very similar systemic velocities, we con- 
clude that these stars are cluster members, but that their veloc- 
ities are off by -39 + 3 km s"' . We suspect that this systematic 
shift may have been introduced by the correction that we ap- 
plied using the night-sky emission line. As this sky Une was 
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Fig. A.2. Normalized spectra of the stars on the slit. The spectra 
are ordered from the star with the lowest V-band magnitude, 
star D, to the brightest, star R. Each consecutive spectrum is 
shifted by one unit. 

located redwards of the reddest wavelength calibration line, 
the polynomial could have introduced this systematic shift. 
However, as mentioned above, we do not expect that this in- 
fluences relative velocities between different observations and 
different stars, since the shape of the polynomial did not vary 
between the different observations. 

The one remaining issue is that of the systemic radial ve- 
locity of star C, which is different from that of the other stars. 
Star C may not be a cluster member, as it does not coincide with 
the cluster main sequence. Furthermore, its spectral features, as 
displayed in Fig. lA.2l are distinct from those of the other stars. 
In this figure, the normalized spectra of the stars are sorted in 
order of increasing V-band magnitude, so that when the stars 
are cluster members their spectra should be ordered on spec- 
tral type and their spectral features should change accordingly. 
Approximate spectral types for these stars were determined by 
comparing the spectra to those in the atlas of Le Borane et al. 

We see that both the spectral features and the spectral 
type of star C show more resemblance with that of star D than 
they do to star R or A. From this we conclude that star C is not 
a cluster member and therefore its radial velocity may differ 
from that of the cluster. 

Appendix B: Is the binary associated with 
NGC 6752? 

We use our observations and the properties of the white dwarf 
and the pulsar that we derived from them, to test if these sup- 
port the assumption that PSR J 19 11 -5958 A is a member of the 
globular cluster NGC 6752. 



B.1. Velocities 

If the pulsar binary is a member of the globular cluster, the 
systemic radial velocity of the binary should be compatible 
with the radial vel ocity of the cluster, -32.0 + 1.6kms ' 
dPubath et allll997l) . 

The absolute systemic velocity of the pulsar binary is y = 
-18 + 6 km s"' . This differs from the radial velocity of the clus- 
ter by 14 + 6kms"' and is consistent with the systemic velocity 
of the binary j^v = 13 ± 6kms ' which we determined from 
the fit of the white dwarf velocities relative to those of star R. 

The velocity difference needs to be corrected for the grav- 
itational redshift of the white dwarf and the reference star R. 
Using the mass and radius of the white dwarf determined in 
§ 14.31 and the mass and radius of star R determined from the 
its absolute magnitude and the models bv lOirardi et al.l (l200(]l) 
plotted in Fig.|2l we obtain 2.7 km s ' and 0.5 km s ' for the 
white dwarf and star R, respectively. This reduces the velocity 
difference to 11 + 6kms"', amounting to about a 2cr differ- 
ence between the velocity of the binary and the velocity of the 
cluster. A part of this velocity difference may be due to the dis- 
persion in the velocity between the pulsar binary and the ref- 
erence star R. From the scatter in the velocities of stars R, A 
and D, we estimate that the velocity dispersion is of the order 
of ~2kms"^ 

Also important is the local escape velocity at the pul- 
sar position relative to the cluster center. To estimate the es- 
cape velocity Vesc = 2GM/r at the projected distance of 
PSR J1911-5958A, we compute the mass M(rj_) inside a 
sphere of radius using relation A3 from the simplified cluster 
model presented in the Appendix of Freire et al. (2003). Here, 
we used a core-radius of = 6'.'7 (Lugger et alill995l) . a dis- 
tance of d = 4.14kpc dGratton et al. . 2003), a ce ntral v eloc- 
ity dispersion of (T^(Q) = 4.5 km s ' (Dubath etail ll997l) and 
rj_ = 637 to obtain M(ri) = 27xlO^M0 andvesc = 5.7km s"'. 
This velocity should be regarded as an upper limit since only 
the projected distance to the cluster center is known and not 
the actual distance = r]_ + r^, with ry the distance along the 
line-of-sight towards NGC 6752. 

We conclude that, taking into account the velocity range 
expected due to the velocity dispersion and the escape ve- 
locity, the systemic velocity of PSR J1911-5958A is only 
marginally consistent (at the 2cr level) with the radial velocity 
of NGC 6752. 



B.2. The radius of the white dwarf 

From § 13 we found that the white-dwarf radius determined 
from the effective temperature and surface gravity is outside 
the Icr uncertainty on the radius derived if the white dwarf is 
assumed to be at the distance of NGC 6752. This suggests that 
the pulsar binary is not a member of the globular cluster. There 
may be additional uncertainties in some of the parameters that 
were used. Here we discuss some of the sources of uncertainty 
that may explain the discrepancy in the white-dwarf properties. 
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B.2.1. Distance, reddening and magnitudes 

As a result of the discrepancy in the white-dwarf radius, there 
is a discrepancy in the distance to the white dwarf. The dis- 
tance modulus (m - M)y = 12.66 ± 0.46 derived from the 
combination of the atmospheric properties of the white dwarf 
and the mass-radius relation (§ I4.4> is only marginally con- 
sistent with the distance modulus {m - M)v = 13.24 ± 0.08 
(Gratton et al., 2003) determined for NGC 6752. Though there 
is a considerable spread in the distance modulus measurements, 
ranging from (m-M)y = 13.17 + 0.13 (c/ = 4.0 ± 0.3 kpc) from 
white-dwarf cooling sequence fitting (Renzini et al., 1996), to 
(m - M)v = 13. 32 + 0.04 (d ^ 4.31 + 0.08 kpc) from main- 
sequence fitting fcarretta et al.l l2000l) . no distance determina- 
tion places NGC 6752 so close. Similarly, the spread in red- 
dening mea surements, Eg^v - 0.035 to 0.046 JGratton etal.L 
I2003[l2005l) . does not help to explain the radius discrepancy. 

It also seems unlikely that our photometry is in error by as 
much as the 0.6 mag or more needed to match the distance and 
hence the white-dwarf radius. Our UBV magnitudes are con- 
sistent with those given by Ferraro et al. (2003), but have sig- 
nificantly smaller uncertainties. The presence of photometric 
B and V standards on the science images removed any uncer- 
tainties on the airmass dependence. Also, the Z?-band calibra- 
tions taken on two different nights were consistent with each 
other, having only a 0.01 mag difference. Finally in § 12.21 we 
found that the white-dwarf companion is not variable (down to 
0.02 mag) and this eliminates the possibility that the UBV pho- 
tometry was taken at an extremum in white-dwarf brightness. 

B.2.2. Line broadening 

As was found in § 14.21 a lower surface gravity, which would 
imply a lower mass and a larger radius, was found when the 
spectral resolution was decreased to 6.5 A. Though this is con- 
siderably larger than the 4.5 A determined from the width of the 
lines in the spectrum of star R, the Balmer lines in the spectrum 
of the white dwarf may be broadened. In order to broaden the 
lines from 4.5 A to 6.5 A a velocity smearing of ~ 300 kms"' 
is required. 

One source of broadening is due to the fact that in a single 
observation, the 2470 s exposure time covers about A(p = 3.5% 
of the 20 h orbit. This introduces a maximum change in velocity 
(at = 0.25 and <f) = 0.75) of about 2nA4>Ky^D ^ 50kms"', 
which is much less than required'. 

Another source of broadening could be due to rotation. To 
estimate the rotational velocity of the white dwarf, we follow 
the reasoning used in van Kerkwiik & Kulkarni ( 1995) to ex- 
plain the variations seen in the spectrum of the white-dwarf 
companion to PSR B0655+64. During the period of mass trans- 
fer, the progenitor of the white dwarf, a giant, filled its Roche 
lobe and tides ensured the system was synchronized and circu- 
larised. When mass transfer ceased and the pulsar companion 
started to contract to a white dwarf, the tides became inefficient 
and the rotational evolution of the companion was determined 

' It also causes a reduction in inferred radial-velocity amplitude by 
a factor sin(;rA0)/;rA(/i = 0.998. This is sufficiently small that we have 
chosen to ignore it. 



by conservation of angular momentum. Thus, the rotational pe- 
riods scale inversely with the moments of inertia. 

The moment of inertia of the progenitor can be separated 
in that of the core and the envelope; /piog - /core + hm with 
/core = kl^rcMcor^Rlow ^nd /env = k^^^M^^^Rl^^, whcrc k is the 
gyration radius. As the progenitor fills its Roche-lobe of radius 
Rl, we have /?env - Rl- After mass transfer, when the envelope 
has contracted to form the white dwarf, the white dwarf has 
/wD = ^WD^WD^WD- Assuming that /core ~ /wD and ignoring 
diff'erences in the gyration radii, conservation of angular mo- 
mentum gives Qrot/Qorb - 1 + M^^^rI/M^^rI^^. Here, two 
assumptions lead to an overestimate of the spin up; the enve- 
lope will be more centrally concentrated than the white dwarf, 
such that fcenv < ^WD, whilc tidal dissipation will be important 
in the initial stages of contraction. On the other hand, the hot 
core of the progenitor will be larger than the white dwarf, so 
/core > /wD (though generally this effect should be small, since 

in most cases /core ^ /env)- 

For a white dwarf with a mass of 0.17-0. 1 8 Mq and a radius 
of 0.042-0.058 /?0, and for the observed mass ratio of ^ = 7.36 
and orbital period of 20 h, the Roche-lobe radius of the pro- 
genitor is about Rl - 0.96 /Jq- After cessation of mass trans- 
fer, the remaining envelope has a mass of about 0.01 Mq. With 
these values we estimate that the rotational period of the white 
dwarf is about 20-30 times faster than the orbital period, so 
fiot ~ 1 to 0.6 h. In that case, the rotational velocity v'rotsin/ 
would be of the order of 50-100 km s ' . This is smaller than the 
300 km s"' estimated above, and since our assumptions likely 
led us to overestimate the rotational velocity, it seems unlikely 
that rotational broadening alone could explain the discrepancy 
between the surface gravity inferred from the spectrum and that 
inferred from the radius assuming that PSR J191 1-5958A is a 
member of NGC 6752. 



6.3. Constraints from radio timing 

b'Amico et 11111200211 give two arguments for the association 
of PSR J191 1-5958A with NGC 6752. The first one is that it 
was discovered in an observation of a globular cluster: the ded- 
icated globular cluster observations with Parkes together cover 
a relatively small area compared to the whole sky, and the de- 
tection of a rare millisecond pulsar in this area suggests mem- 
bership of the cluster. The problem is that the dedicated clus- 
ter observations are much deeper than most other pulsar ob- 
servations, and that the number of millisecond pulsars at the 
flux level of PSR J1911-5958A, 0.22 mJy at 1400 MHz, and 
their distribution on the sky, are not known. Estimates based 
on extrapolation are uncertain. The accuracy of the estimate by 
D' Amico et al. of a 10'^ probability for a chance coincidence is 
therefore not clear. lD' Amico et alT ( l2002h used the center beam 
of the Parkes multibeam receiver in their disco very observa- 
tion. The diameter of that beam is about 14' ( Manchester et al.L 
iool), and thus any pulsar detected in the globular cluster sur- 
vey must lie within 7' from the cluster center. We conclude that 
the argument from the probability of chance coincidence is less 
solid than the number 10"^ might suggest. 
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The second argument of iD'Amico et alJ ( l2002i) is that 
PSR J1911-5958A has a dispersion measure DM = 
33.7 pc cm"^^, almost the same as the dispersion measure DM = 
33.3 pccm"^ of the three pulsars in the cluster core, two of 
which certainly belong to the cluster as proven by a negative 
period derivative. According to the iTavlor & CordesI lll993h 
model, the maximum dispersion measure in the direction of 
NGC6752 and PSR J1911-5958A is DM,™x ~ 42pccm-^ 
Since the DM values of the pulsars in the core of NGC 6752 
(and of PSR J1911-5958A) ai-e less than this, the pulsar 
would be almost at the distance of the cluster. However, the 
maximum to the dispersion measure arises because the elec- 
tron layer in our Galaxy has a finite scale height of 0.5 to 
1 kpc, and this implies that all pulsars above the electron 
layer have the same dispersion measure in the same direction 
(Bhattacharva& Verbunt, 1991). Since NGC 6752 is at a dis- 
tance of d - 4. 14 kpc and a Galactic latitude of b - -25°6, 
it is well above the electron layer, and its observed DM must 
be equated with the actual DM^ax in that direction. (This is 
compatible with the uncertainty of about V4DM^^ in the 
model value of DM,„ax, which gives 13pccm"^ for DM^ax ~ 
42pccm"-'; Nelemans et al. 1997). We conclude that the dis- 
persion measure of PSR J 191 1-5958 A does not prove that it is 
a member of the cluster, but merely that it is above the galactic 
electron layer, i.e., at a distance d ^ 2.4 kpc. 

B.4. White dwarf models 

Finally, we cannot discard the possibility that the white-dwarf 
models themselves are uncertain. This can already be seen from 
the mass-radius relations shown in Fig.|6l These relations are 
for the observed temperature of Teff = 10090 K, but, for a given 
mass of say, 0.20 Mq, the predicted radii show a spread of about 
0.01 Ro. 

A part of this problem is the lack of low-mass, helium-core 
white dwarfs for which accurate parameters have been deter- 
mined and which could be used to calibrate the evolution, cool- 
ing and atmospheric properties of these systems. 

B.5. Summary and prospects 

Summarizing the results of this section, we see that, at face 
value, the systematic radial velocity and radius estimate indi- 
cates that PSR J1911-5958A is not a member of NGC6752. 
Furthermore, we have argued that the similarity in dispersion 
measure for all all five pulsars located towards NGC 6752 does 
not necessarily imply that they are all at the same distance. 

However, our observations do not conclusively rule out the 
membership of the pulsar and the globular cluster either. All 
constraints are consistent at the 2cr level, and the inconsisten- 
cies of the constraints may be removed when we take into ac- 
count that there is an allowed range in radial velocity due to the 
velocity dispersion, that there may be broadening of the B aimer 
lines in the spectrum of the white dwarf and that there are pos- 
sible uncertainties in the white-dwarf models themselves. As 
such, it is for future observations to decide between either pos- 
sibility. 



